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We show that given a nonvanishing particular solution of the equation 



the corresponding differential operator can be factorized into a product 
of two first order operators. The factorization allows us to reduce the 
equation (1) to a first order equation which in a two-dimensional case is the 
Vekua equation of a special form. Under quite general conditions on the 
coefficients p and q we obtain an algorithm which allows us to construct in 
explicit form the positive formal powers (solutions of the Vekua equation 
generalizing the usual powers (2 — 2;o)",n = 0,l,...). This result means 
that under quite general conditions one can construct an infinite system 
of exact solutions of (1) explicitly, and moreover, at least when p and q are 
real valued this system will be complete in ker(divpgrad +g) in the sense 
that any solution of (1) in a simply connected domain n can be represented 
as an infinite series of obtained exact solutions which converges uniformly 
on any compact subset of SI. 

Finally we give a similar factorization of the operator (divpgrad +g) in 
a multidimensional case and obtain a natural generalization of the Vekua 
equation which is related to second order operators in a similar way as its 
two-dimensional prototype does. 



1 Introduction 

Consider the one-dimensional stationary Schrodinger equation 



Abstract 



(divpgrad +g)it = 0, 



(1) 




(2) 
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It is well known that given a nonvanishing particular solution of ||2Jl, the 
Schrodinger operator can be factorized 




(3) 



and as a consequence the general solution of (0) can be obtained. 

In )18j it was shown that a similar factorization of a stationary Schrodinger 
operator is available in two dimensions. One of the main results of the present 
work is the factorization of a more general two-dimensional elliptic operator 
(divpgrad+g) in a form similar to Q. The factorization allows us to reduce 
the equation 



to a Vekua equation of a special form. This Vekua equation becomes bicomplex if 
p or q are complex functions. One complex component of a solution of the Vekua 
equation (its real part when p and q are real valued) necessarily satisfies equation 
(0J and the other satisfies an associated equation having the form of 10} but with 
different coefficients p and q. This situation is a generalization of the fact that 
real and imaginary parts of an analytic function are harmonic, and likewise for 
any harmonic function in a simply connected domain its harmonic conjugate 
can be constructed, we obtain explicit formulas for constructing "conjugate" 
solutions of the associated equations of the form In the case p = I and 
q = these formulas turn into the well known from complex analysis formulas 
for construction of conjugate harmonic functions. 

In JH| we established that under quite general conditions the positive formal 
powers corresponding to the Vekua equation (the pseudoanalytic functions gen- 
eralizing the usual powers (z — zo)", n = 0, 1, . . .) can be constructed explicitly. 
Here we develop this result and obtain that by a given nonvanishing particular 
solution of equation 10} under quite general conditions one can construct an 
infinite system of exact solutions of (@} explicitly, and moreover, at least when 
p and q are real valued this system will be complete in ker(divpgrad-l-g) in the 
sense that any solution of Q in a simply connected domain f2 can be represented 
as an infinite series of obtained exact solutions which converges uniformly on 
any compact subset of il. 

In the final part of the present work we obtain a factorization of the operator 
(divpgrad-|-(?) in a multidimensional situation and reduce 10} to a first order 
equation which generalizes the Vekua equation. 

2 Preliminaries 

A bicomplex number has the form 



where Qi = 90 + iqii Q2 = 92 + iqs, qj G M, j = 0, 3; — k"^ = —1, ik — ki. 



(divpgrad-|-(7)u = 



(4) 



9 = Qi + Q2k 
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We will say that Qi and Q2 are complex components of the bicomplex num- 
ber q. Denote q = Qi — Q2k. The corresponding conjugation operator we denote 
by C: Cq — q. We will say that q is scalar if Q2 = 0. That is in this case q is 
a usual complex number. Q2 is called the vector part of Q. Sometimes we use 
the notation: Qi = Sc{Q), Q2 = Vec{Q). 

The set of bicomplex zero divisors, that is of nonzero elements q — Qi + Q2k-, 
{Qi, Q2} C C such that 

qq = (Qi + Q^k) (Qi - Qafc) = (5) 

we denote by S. 

We will consider the variable z = x + ky, where x and y are real variables 
and the corresponding differential operators 

(h=^{dx + kdy) and = ]^{dx- kdy) 

Notation W-g or Wz means the application of or dz respectively to a bicomplex 
function W{z). 

Note that if we consider the equation 

dz^ = $ (6) 

in a whole complex plane or in a convex domain, where $ = (f>i + fc$2 is a 
given bicomplex valued function such that its scalar part $1 and vector part $2 
satisfy the equation 

ay$i + 9,$2 = 0, (7) 

then there exists a scalar solution of (jSJ which can be reconstructed up to an 
arbitrary scalar constant c in the following way 



(8) 



\J xa J yo 

where (xq, yo) is an arbitrary fixed point in the domain of interest. 
By A we denote the integral operator in 

Amx,y)^2(f <i>i(77,y)d77- T <i>2(a;o, O^^f 

\J xo J yo 

Note that formula ijHJl can be easily extended to any simply connected domain 
by considering the integral along an arbitrary rectifiable curve T leading from 
{xo,yo) to {x,y) 

ip{x, y)^2(^J ^idx - $2dy^ + c. 

Thus if $ satisfies 0, there exists a family of scalar functions (p such that 
dzf = given by the formula (p = A[^]. 

In a similar way we define the operator A corresponding to dz- 

Amx,y)^2( f <^>i{7j,y)dv+ f ^2{xo,0d^ 

\Jxo Jyo 
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3 Solutions of second order elliptic equations as 
scalar parts of bicomplex pseudoanalytic func- 
tions 

Consider the equation 

(_A + I.) / = (9) 

in some domain C R^, where A — + v and / are complex valued 
(in our terms scalar) functions. We assume that / is a twice continuously 
differentiable function. 

Theorem 1 Let j he a nonvanishing in fl particular solution of Then for 
any complex valued (scalar) function ip G C^(f2) the following equalities hold 

\{A-iy)ip= (^+ jc) {dz - jc) ^ - + (^(h- V- 

(10) 

Proof. Consider 

(3.+ kc) (a - Ifc) ^ . li, - - * (^) , 

= -(A^--^^) = -(A-z.)^. 
Apphcation of C gives us the last part of |[TU|l . ■ 

In the case of a real valued potential v this theorem was proved in |18| . 
The following statement is known in a form of a substitution (see, e.g., |23)'). 
Here we formulate it as an operator relation. 

Proposition 2 Let p and q be complex valued functions, p G C'^{fl) and p 
in fl. Then 

divpgrad+g =p^/^(A - r)p^/^ in il, (11) 

where 

_ ApV^ _ q 

pl/2 p 

Proof. The easily verified relation 

divpgrad = pi/2(A-^^)pi/2 (12) 

is well known (see, e.g., |22)- Adding to both sides of 1)12(1 the term q (and 
representing it on the right-hand side as p^/^ {q/p)p^^'^) gives us (fTT|l . ■ 
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Theorem 3 Let uq be a nonvanishing in Q particular solution of the equation 

(diY pgTa,d+q)u = in fl. (13) 

Then under the conditions of proposition\^for any complex valued (scalar) con- 
tinuously twice differ entiable function the following equality holds 

i(divpgrad+g)^=pi/2(^a, + ^C^ (^ck ^ p^^^, (14) 

where 

f^p'^W (15) 

Proof. This is based on (|10|l . From (|ll|l we have that if uq is a solution of p3|l 
then the function (|15|l is a solution of the equation 

(A - r)f ^ 0. (16) 

Then combining ((TT1) and ifTUI) we obtain ifTH) . ■ 

Remark 4 According to iTT^) . A — r = /^^ div grad/^^ w/iere / is a solution 
of nib]} . Then from we have 

divpgrad+q =//V"^div/2grad/"V/2. (17) 

Taking into account il we obtain 

div p grad +<z = Wg^ ^ div puQ grad ^ m . 

Remark 5 Let q = 0. Then uq can be chosen as uq = 1. Hence fj^l ) gives us 
the equality 

idiv(pgrad^) (^5^ + ^^^) b^^V). 

In what follows we suppose that in il there exists a nonvanishing particular 
solution of p3|l which we denote by uo- 

Let / be a scalar function of a; and y. Consider the bicomplex Vekua equation 

= jW in O. (18) 

Denote Wi = ScW and W2 = Vec W. 
Remark 6 111^ Equation can be written as follows 

fikif-'Wi) + kf-'ck{fW2) = 0. (19) 
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Theorem 7 Let W = Wi + be a solution of U^. Then U = f~'^Wi is a 
solution of the equation 

div(/2vC/) =0 inn, (20) 
and V — fW2 is a solution of the equation 

divir^VV) = inn, (21) 
the function Wi is a solution of the stationary Schrddinger equation 

-APFi +riVFi = inn (22) 
with Ti = Af / f, and W2 is a solution of the associated Schrddinger equation 

-AW2 + raVFa = inn (23) 
where r2 - 2{V ff / f^ - n and {V ff = fl + f^ 

Proof. To prove the first part of the theorem we use the form of equation H18|l 
given in Remark El Multiplying H19|) by / and applying gives 

dJ^ch{f-^Wi) + -^A{fW2) =0 

from where we have that So [dzf^dz (/~^Wi)) = which is equivalent to H20|l 
where U = f'^Wi. 

Multiplying ifTi^ by /^^ and applying dz gives 

i A (/-^W^i) + kdzf-^^ {fW2) - 

from where we have that Sc (Sz/^^c^ (/W2)) = which is equivalent to ^(T\\ 
where V = fW2- 

From H12|l we have 

(A - n) Wi = /~i div(/2v {f-^Wi)). 

Hence from the just proven equation (|20|l we obtain that Wi is a solution of 
(EH). 

In order to obtain equation H23I) for W2 it should be noted that 
/div(r2v(/M^2)) = (A-r2)W2. 

■ 

In the case of a real valued function / the relation between solutions of H18|l 
and equations 12UII . (|21|) was observed in |20|, and between solutions of (|18|) and 
equations HSJ in llBj- 
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Remark 8 Observe that the pair of functions 

F = f and G^j (24) 

is a generating pair for This allows us to rewrite ^18(1 in the form of an 

equation for pseudoanalytic functions of second kind 

Vzf + = 0, (25) 

where and ip are scalar functions. If ip and %p satisfy 1^25}) then W — (pf + '4'j 
is a solution of and vice versa. 

Denote w — ip + ijjk. Then from \25}) we have 

{w + w)-f + {w ~ "^hj = 0' 
which is equivalent to the equation 

1 - /2_ 

^~ = TT72 ^~ ^'^^^ 

The relation between \26\) and \20\) . I^21\) in the case of a real valued function f^ 
was observed in "^Z and resulted to be essential for solving the Calderon problem 
in the plane. 

Theorem 9 Let W = Wi + W2k be a solution of fl^) . Assume that f = p^/^uo, 
where uq is a nonvanishing solution of f J/)]) in Vl. Then u = p~^^^Wi is a 
solution of II 1,'^} in Q, and v — p^/'^W2 is a solution of the equation 

(div - grad +qi)v = inil, (27) 
P 



where 

p \p \ p uo / \ uo J 



Proof. According to theorem[7| the function / ^Wi is a solution of (|20|) . From 
p7|l we have that 

(divpgrad+g) (p-^/^Wi) = f-^ div{f^W {f-^Wi)) 

from which we obtain that u = p^^^^Wi is a solution of 

In order to obtain the second assertion of the theorem, let us show that 

pi/2(divigrad+qi)(pi/V) - /div(/-2v(/^)) 
P 
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for any scalar £ C'^{VL). According to p2|l . 

/div(/-2v(/^)) = (^A - V={^- r2) V. 

Straightforward calculation gives us the following equality 

A/'^ _ 3 /V^y 1 Ap ^ /Vp Vuo\ Auo^^fVuoy 

From the condition that uq is a solution of we obtain the equality 

Auo <? / Vwq \ 



"0 P \ P uo 
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Thus, 



Notice that 

2 



A/-1 _ 3 /^VpV 1 Ap , ^/Vp Vuo\ , 9 , ^ /^Vuo 



Ap-i/2 3 /VpV 1 Ap 



/Vp 








\ p ' 


Uo 1 


P 





p-i/2 4 V p / 2 p 

Then 

Now taking qi in the form (|28|l we obtain the result from ■ 

Theorem 10 JT^' Lei VKi be a solution of ^ in a simply connected domain 
Q. Then the function W2, solution of l^2,S\) such that W — Wi + W2k is a 
solution of IjiSj) . is constructed according to the formula 

W2 = f-'A{kf^ck{f-^Wi)). (29) 

It is unique up to an additive term cf^^ where c is an arbitrary complex constant. 

Given a solution W2 of the corresponding solution W\ of ^2^) such 

that W = Wi + W2A; is a solution of is constructed as follows 

Wi = -fA{kf-^ck{fW2)) (30) 

up to an additive term cf. 

Remark 11 When 1^ = and /o = 1, equalities ^2y\) and HclUji) turn into the 
well known in complex analysis formulas for constructing conjugate harmonic 
functions. 
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Corollary 12 Let U be a solution of \2U\) . Then a solution V of such that 

W = fU + kf-^V 
is a solution of is constructed according to the formula 

V - A{kf^U-). 

It is unique up to an additive complex constant. Conversely, given a solution V 
of \21}) . the corresponding solution U of l^2(A) can be constructed as follows: 

U = -A{kf-^V^). 

It is unique up to an additive complex constant. 

Proof. Consists in substitution of Wi = fU and of W2 — f^^V into and 
Ipp . ■ 

Corollary 13 Let f — p^^^uq, where is a nonvanishing solution of !^13\) in 
a simply connected domain Q, and u be a solution of Then a solution v 

of I27[ l with qi defined by h2Sp such that W = p^^^u + kp^^/'^v is a solution of 
is constructed according to the formula 

V = UQ^A{kpuQdziuQ^u)). 

It is unique up to an additive term cUq^ where c is an arbitrary complex constant. 

Let V be a solution of \2"]\j , then the corresponding solution u of 1) such 
that W — p^^^u + kp^^/'^v is a solution of \lt^) . is constructed according to the 
formula 

u = —uoA{kp^^UQ^dz{uo'^))- 

Proof. Consists in substitution of / = p^^^uo, Wi = p^^^u and W2 = p^^^'^v 
into ^ and m 



4 Some definitions and results from pseudoana- 
lytic theory for bicomplex functions 

4.1 Generating pair, derivative and antiderivative 

Following [5j we introduce the notion of a bicomplex generating pair. 

Definition 14 A pair of bicomplex functions F — Fi +F2k and G = Gi +G2k, 
possessing in Q partial derivatives with respect to the real variables x and y is 
said to be a generating pair if it satisfies the inequality 

Vec(FG) ^0 infl. 

The following expressions are called characteristic coefficients of the pair [F, G) 

FG— — F—G FG— — F-^G 

""^^'^^ ^ " FG-FG ' ^(^'^^ ^ FG^FG ' 
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FGz — F^G FGz — F^G 

A{F,G) = — — =7r' ^{F,G) = 



FG-FG ' ' ' FG-FG 

Every bicomplex function W defined in a subdomain of admits the unique 
representation W = 4'F + ipG where the functions cf) and ip are scalar. 

The {F, G')-derivative W = '^'''^fj ^ of a function W exists and has the form 
W = cPzF + i,zG = Wz- A^F^G)W - B(F^G)W (31) 

if and only if 

(hF + 4>zG = 0. (32) 
This last equation can be rewritten in the following form 

Wz = a^^p^G)W + b(F^G)W (33) 

which we call the bicomplex Vekua equation. Solutions of this equation are 
called {F, G)-pseudoanalytic functions. 

Remark 15 The functions F and G are {F, G)-pseudoanalytic, and F = G = 
0. 

Definition 16 Let {F,G) and {Fi,Gi) be two generating pairs in fi. (Fi,Gi) 
is called successor of (F, G) and {F, G) is called predecessor of (Fi, Gi) if 

«(Fi,Gi) = 0-{F,G) '^'"'d ^(Fi,Gi) = -B(F,G}- 

The importance of this definition becomes obvious from the following state- 
ment. 

Theorem 17 Let W he an {F,G)-pseudoanalytic function and let (Fi,Gi) be 

a successor of {F, G). Then W is an {Fi,Gi)-pseudoanalytic function. 

Definition 18 Let (F, G) be a generating pair. Its adjoint generating pair 
(F, G)* = (F*,G*) is defined by the formulas 

2F 2G 
F* = = — ^— , G = 



FG-FG FG-FG 
The (F, G)-integral is defined as follows 

J Wd(^F,G)Z = F{zi)Sc J G*Wdz + Gizi)Sc j F*Wdz 

where F is a rectifiable curve leading from zq to Zi . 

If = 0F + ipG is an (F, G)-pseudoanalytic function where ^ and V are 
scalar functions then 



r 

J Zn 



Wd^F^G)Z = W{z) - <f>{zo)F{z) - i^{zo)G{z), (34) 



and as F =G = 0, this integral is path-independent and represents the (F, G)- 
antiderivative of W. 
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4.2 Generating sequences and Taylor series in formal pow- 
ers 

Following '5' we introduce the following definitions and results. 

Definition 19 A sequence of generating pairs {(F,„, Gm)}, m = 0, ±1, ±2, . . . 
is called a generating sequence if {Fm+i,Gm+i) is a successor of {Fm,Gm)- If 
(Fo,Go) = (-^iG), we say that {F,G) is embedded in {{Fm,Gm)}- 

Theorem 20 Let {F, G) he a generating pair in 57. Let fii he a hounded do- 
main, f2i C ri. Then {F,G) can he embedded in a generating sequence in 

Definition 21 A generating sequence {{Fm, Gm)} is said to have period /i > 
if {Fm+fi, Gm+n) is equivalent to {Fm,Gm)i that is their characteristic coeffi- 
cients coincide. 

Let W be an {F, G)-pseudoanalytic function. Using a generating sequence 
in which {F, G) is embedded we can define the higher derivatives of W by the 
recursion formula 

dz 

Definition 22 The formal power Zm^ct, zq; z) with center at zq G O,, coefficient 
a and exponent is defined as the linear combination of the generators Fm, Gm 
with complex constant coefficients \, fi chosen so that XFm{zo) + fJ-Gmizo) — o- 
The formal powers with exponents n = 1,2,... are defined by the recursion 
formula 

Z(;'+i)(a,zo;^) = (n+l) rzi"ji(a,zo;CKF™,G„)C- (35) 

J Zo 

This definition implies the following properties. 

1. Zm \a, zq] z) is an {Fm, G'm)-pseudoanalytic function of z. 

2. If a' and a" are scalar constants, then 

Zl^^a' + ka", Zo; z) = a'Z^)(l, zq; z) + a" Z^^^k, zo; z). 

3. The formal powers satisfy the differential relations 

d(Fr,,,G^)Zm\a,Zo;z) ^(n-l)/ N 

^ = nZ'm+i[a, Zo; z). 

4. The asymptotic formulas 

Zi")(a,zo;z)^a(z-zo)", z ^ zq (36) 

hold. 
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Assume now that 

oo 

W{z) = Y,Z^^Ha,z^;z) (37) 

where the absence of the subindex m means that ah the formal powers corre- 
spond to the same generating pair {F,G), and the series converges uniformly 
in some neighborhood of zq. It can be shown that the uniform limit of pseu- 
doanalytic functions is pseudoanalytic, and that a uniformly convergent series 
of {F, G)-pseudoanalytic functions can be {F, G)-differentiated term by term. 
Hence the function W in (|37l) is {F, G)-pseudoanalytic and its rth derivative 
admits the expansion 

oo 

W^"\z) = ^ n(n - 1) • • • (n - r + l)Z("-'-)(a„, zq; z)- 



From this the Taylor formulas for the coefficients are obtained 

W^["l(zo) 



(38) 



Definition 23 Let W{z) be a given (F, G) -pseudoanalytic function defined for 
small values of \z — zq\. The series 

oo 

5]z(")(a,zo;z) (39) 

with the coefficients given by is called the Taylor series ofW at zq, formed 
with formal powers. 

The Taylor series always represents the function asymptotically: 

N 

W-(z)-^Z(")(a,Zo;^) = 0(1^-^01"^+'), z^zo, (40) 

ri=0 

for all N. 

If the series (|39f) converges uniformly in a neighborhood of zo, it converges 
to the function W. 



4.3 Convergence theorems 

The statements given in this subsection were obtained by L. Bers and S. Agmon 
and L. Bers. Their proof in a usual complex case was based on the so called sim- 
ilarity principle. The similarity principle in general is not valid in a bicomplex 
situation. Here we correct the corresponding statement which unfortunately in 
1^ was formulated with a mistake. 
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Theorem 24 (Similarity principle) Let w be a regular solution of in a 
domain fl such that its values are not zero divisors at any point. Then the 
bicomplex function ^ — w ■ e'^, where 



is a solution of the equation 9^$ = in Q. 

The proof is completely analogous to that for a complex case (see |25|^. 

Up to now this is an open question how one can guarantee that a solution of 
the bicomplex Vekua equation (|33|l be different from a zero divisor at any point. 
This is why the proof of the following statements is valid in a usual complex 
situation, and in a bicomplex case their validity should be investigated. 

In the case when the coefhcients in (|33|l are usual complex functions (with 
respect to k) the following theorems regarding the convergence of formal Taylor 
expansions are valid. 

Theorem 25 The formal Taylor expansion 1^3 9\) of a pseudoanalytic func- 
tion in formal powers defined by a periodic generating sequence converges in 
some neighborhood of the center. 

Definition 2G \Bj A generating pair {F, G) is called complete if these functions 
are defined and satisfy the Holder condition for all finite values of z, the limits 
F{oo), G{oo) exist, Vec{F{oo)G{ooj) > 0, and the functions F{l/z), G{l/z) 
also satisfy the Holder condition. A complete generating pair is called normal- 
ized if F{oo) — I, G(oo) = k. 

A generating pair equivalent to a complete one is complete, and every com- 
plete generating pair is equivalent to a uniquely determined normalized pair. 
The adjoint of a complete (normalized) generating pair is complete (normal- 
ized). 

From now on we assume that (F, G) is a complete normalized generating pair. 
Then much more can be said on the series of corresponding formal powers. We 
limit ourselves to the following completeness results (the expansion theorem and 
Runge's approximation theorem for pseudoanalytic functions). 

Following [5] we shall say that a sequence of functions Wn converges normally 
in a domain Q if it converges uniformly on every bounded closed subdomain of 



Theorem 27 Let W be an (F, G) -pseudoanalytic function defined for \z — zo\ < 
R. Then it admits a unique expansion of the form W{z) = X]^o ^''"H'^nj ^oj -z) 
which converges normally for |z — zo| < 9R, where 9 is a positive constant 
depending on the generating sequence. 




a{F,G){z) + b[F.G){z)^^ ifw{z)^0, z £ fl, 
a{F,G){z) + b{F,G){z) ifw{z)^0, zE^ 
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The first version of tliis theorem was proved in PP. We follow here [7]. 

Remark 28 Necessary and sujficient conditions for the relation 9^1 are, 
unfortunately, not known. However, in the following sufficient conditions 
for the case when the generators {F, G) possess partial derivatives are given. 
One such condition reads: 

I / M I ^ / M ^ Const 
\F^{z)\ + \G^{z)\ < 

1 + \z\ 

for some s > 0. Another condition is 

f f (\F-\'-' + \F-f+' + \G-\'-' + \G-\'+')dxdy<^ 

J J\z\<oo ^ ' 

for some < e < 1. 

Theorem 29 ^ A pseudoanalytic function defined in a simply connected do- 
main can be expanded into a normally convergent series of formal polynomials 
(linear combinations of formal powers with positive exponents). 

Remark 30 This theorem admits a direct generalization onto the case of a 
multiply connected domain (see ^). 

In posterior works [221 ) deep results on interpolation and on the 
degree of approximation by pseudopolynomials were obtained. For example, 



Theorem 31 Let W be a pseudoanalytic function in a domain il bounded 
by a Jordan curve and satisfy the Holder condition on dfl with the exponent a 
(0 < a < 1). Then for any e > and any natural n there exists a pseudopoly- 
nomial of order n satisfying the inequality 

. M Const — 
\W{z) - P„(z)| < for any z en 

where the constant does not depend on n, but only on e. 
The primary aim of the next section is to show that: 

1. all the mentioned results are of immediate application to the equation 

2. in many practically important situations the generating sequence and con- 
sequently the formal powers n = 0, 1, . . . can be constructed explic- 
itly. 
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5 Complete systems of solutions for second or- 
der equations 

In what follows let us suppose that the scalar function / is defined in a somewhat 
bigger domain with a sufficiently smooth boundary. Then we change the 
function / for z S rig\f2 and continue it over the whole plane in such a way that 
/ = 1 for large \z\ (see [7 ). In this way the generating pair {F,G) — {f,k/f) 
becomes complete and normalized. 

Then the following statements are direct corollaries of relations established 
in section O between pseudoanalytic functions (solutions of (|18|l ) and solutions 
of second order elliptic equations, and convergence theorems from the previous 
section. 

Definition 32 Let u{z) be a given solution of the equation Hlci\} defined for 
small values of \z — zq\, and let W{z) be a solution of constructed according 
to corollarv MSl such that ScW — p^/'^u. The series 

oo 

p-i/2(z)^ScZ(")(a„,zo;^) 

n=0 

with the coefficients given by is called the Taylor series of u at zq, formed 
with formal powers. 

In the rest of this section we assume that all the coefficients in second order 
equations considered in section 13 are real valued functions and the particular 
nonvanishing solution uq of (|13() is real valued as well. 

Theorem 33 Let u{z) be a solution of defined for \z — zq\ < R. Then it 
admits a unique expansion of the form 

oc 

7/(z)=p-i/2(z)^ScZ(")(a„,zo;2) 
which converges normally for |z — zo| < i?. 

Proof. This is a direct consequence of theorem 12 71 and remark ESI Both neces- 
sary conditions in remark [2H1 are fulfilled for the generating pair H24|) . ■ 

Theorem 34 An arbitrary solution of defined in a simply connected do- 
main where there exists a nonvanishing particular solution uo can be expanded 
into a normally convergent series of formal polynomials multiplied byp"^/"^. 

Proof. This is a direct corollary of theorem [5^1 ■ 

More precisely the last theorem has the following meaning. Due to Property 
2 of formal powers we have that Z'"^ (a, zg; z) for any Taylor coefficient a can be 
easily expressed through Z("^(l,zo;z) and Z'^"'\k^ zq; z). Then due to theorem 
EHlany solution W of H18(l can be expanded into a normally convergent series 
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of linear combinations of Z^'^\l, zq; z) and Z^") (fc, zo; z). Consequently, any 
solution of (|13|l can be expanded into a normally convergent series of linear 
combinations of scalar parts of Z^^\l^ Zq] z) and Z^'^\k,zo]z) multiplied by 

Obviously, for solutions of H13|l the results on the interpolation and on the 
degree of approximation like, e.g., theorem]^ are also valid. 

Let us stress that theorem 1341 gives us the following result. The functions 

{p-^'\z)ScZ^-^\l,zo-z), p-i/2(z)ScZ(")(fc,zo;z)r (41) 

represent a complete system of solutions of (|13f) in the sense that any solution 
of H13() can be approximated arbitrarily closely by a normally convergent series 
formed by functions (|41(l in any simply connected domain Q, where a positive 
solution of (|13l) exists. Moreover, as we show in the next section, in many 
practically interesting situations these functions can be constructed explicitly. 



6 Explicit construction of positive formal pow- 
ers 

The book (see also |10| Supplement to Chapter 4]) contains explicit formulas 
for calculation of positive formal powers in the case when F and G have the 
form 

In ^H] the class of generating pairs for which the generating sequence and 
hence the corresponding formal powers can be constructed explicitly was sub- 
stantially extended. For the generating pair of the form (|24ll it is possible when 
/ fulfils the following condition. 



Condition 35 (Condition S) '18t Let f be a scalar function of some real vari- 
able p : f ^ f{p) 

it by s{p) = 



able p : f ^ f(p) such that the expression is a function of p. We denote 



Besides the obvious example of any harmonic function p and as a consequence 
of p being a Cartesian variable ov p — argz — arctan(?//x), there are many 
other practically important examples of p satisfying Condition S. An important 
example is p{x, y) = \J x^ + y^. In this case s(p) — ^. The parabolic coordinate 

p(x,y) — y/x^ + y^ + X also fulfills Condition S: s{p) — EUiptic coordinates 
fulfil Condition S as well (see |18p. 

Denote by S an antiderivative of s with respect to p. 

Theorem 36 ]18f Let f be a scalar function of a real variable p satisfying 
Condition S and let the function (p — ke^^ pz have no zeros and be bounded 
in Q. Then the generating pair {F,G) with F = f and G = k/f is embedded 
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in the generating sequence {Fm,Gm), ™ = 0, ±1,±2,... with F — (p"^F and 
G = ip'^G. 

This result opens the way for expUcit construction of positive formal powers 
for the equation (|18|l and as a consequence of the complete system of solutions 
ijIT)) for equation ((T^ . 

Some examples of explicitly constructed formal powers were given in jl8) . 
Here we show another quite simple but illustrative example. 

Example 37 Consider the Helmholtz equation 

(-A + c^)u = (42) 

with c being a real constant. Take the following particular solution of 

/ = e'^'^ . Let us construct the first few corresponding formal powers with center 

at the origin. We have 

Z(")(l,0;z) = e^y, Z'^°\k,0;z) = ke'^y , 
Z«(l,0;z) = xe- + ^(i)(fc,0; z) = -'^^ + fcxe"-, 



Z(2)(l,0;z)= (a;2-^)e^^ + 



sinh(cy) 2kx sinh(c?;) 



0; z) = + k (^[x' + ^) e-^y - , . . . . 

It is a simple exercise to verify that indeed the asymptotic formulas (|36(l hold. 
Now taking scalar parts of the formal powers we obtain a complete system of 
solutions of the Helmholtz equation: 

f \ cy ( \ cy / N sinh(cy) 

ui{x,y) = ey, U2(x,y)^xey, U3{x,y) = , 

c 

/ 2 y\ cy , sinh(c?;) 2a;sinh(cy) 

U4{x,y)=[x -- e^H ^ ' uz{x,y) ^ 

V c/ c c 

Formal powers of higher order can be constructed explicitly using a computer 

system of symbolic calculation. For this particular example (together with Maria 

Rosalia Tenorio) Matlab 6.5 allowed us to obtain analytic expressions for the 

formal powers up to the order ten, that gave us the first twenty one functions 

iti, . . . ,M2i. We used them for a numerical solution of the Dirichlet problem 

for the Helmholtz equation with very satisfactory results. For example, in the 

case when f2 is a unit disk with centre at the origin, c — 1 and u on the 

boundary is equal to (this test exact solution gave us the worst precision 

because of its obvious "disparateness" from functions ui,U2 ■ . ■) the maximal 

error max^gn |w(z) — u{z)\ where u is the exact solution and u — J2n=i o-nUn, 

the real constants a„ being found by the collocation method, was of order 10^^. 

A very fast convergence of the method was observed. 
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Although the numerical method based on the usage of explicitly or numeri- 
cally constructed pseudoanalytic formal powers still needs a much more detailed 
analysis these first results show us that it is quite possible that in due time and 
with a further development of symbolic calculation systems it can rank high 
among other numerical approaches, especially for solving equations lO or H13|l 
with rapidly varying coefficients, when finite-difference methods fail. 

7 Reduction of the multidimensional second or- 
der equation to a first order equation 

Here we consider the case of dimension n = 3 and in the final part of this 
section we show that a simple generalization gives us the same results in higher 
dimensions. 

We will consider the algebra ]HI(C) of complex quaternions or biquaternions 
which have the form Q = Qo+ Qi'i + Q2i + Qsk, where {Qk} C C, and i, j, k 
are the quaternionic imaginary units. 

The vectorial representation of a complex quaternion will be used. Namely, 
each complex quaternion Q is a sum of a scalar Qo and of a vector Q: 

Q = Sc(Q) -f Vec(Q) = Qo + Q, 

where Q = Qii + Q23 + Qsk. The operator of quaternionic conjugation we 
denote by Ch- Q — CrQ = Qo ^ Q- We conserve the bar for the quaternionic 
conjugation which should not provoke any confusion with the same notation 
for the conjugation in the first part of the paper because essentially it can be 
considered as the same operation if the bicomplex numbers are considered being 
embedded in IHI(C) in a natural way. 

The purely vectorial complex quaternions (Sc(Q) = 0) are identified with 
vectors from C"^. Note that = ^ < Q, Q > where < •, • > denotes the usual 
scalar product. 

By we denote the operator of multiplication by a complex quaternion 
P from the right-hand side: M^Q = Q ■ P. More information on the structure 
of the algebra of complex quaternions can be found for example in or |21) . 

Let Q be a complex quaternion valued differentiable function of x = (a;, y, z). 
Denote 

DQ = i—Q+i—Q + k—Q. 

ox oy oz 

This expression can be rewritten in a vector form as follows 

DQ = - div Q grad Qo + rot Q. 

That is, Sc(£'Q) = — div Q and Vec{DQ) — grad Qq -I- rot Q. Let us notice that 
= —A. If Qo is a scalar function then DQq coincides with gradQo. 
The following generalization of Leibniz's rule can be proved by a direct 

calculation (see p. 24]). 
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Theorem 38 (Generalized Leibniz rule) Let {P,Q} C C^{G;M{C)), where G 
is some domain in M.^ . Then 

D[P-Q]=D[P]-Q + P-D[Q] + 2{Sc{PD))[Q], (43) 

where 

3 

{SciPD))[Q] :=-^P,a,Q. 

We will actively use the following 

Remark 39 // in r/ieorem[13 Vec(P) = 0, that is P = Pq; then 

D[Po ■ Q] = D[Po] ■ Q + Po ■ D[Q]. (44) 

From this equality we obtain that the operator D + can be factorized as 

follows 

+ i^^)Q ^ P„-'D{Po ■ Q). (45) 

Let G be a complex valued vector such that rotG = 0. Then the complex 
valued scalar function ip is said to be its antigradient if grad(y3 = G. We will 
write if = ^[G]. The operator ^ is a simple generalization of the usual an- 
tiderivative and of the operator A (see section 2), and it defines the function ip 
up to an arbitrary constant. Its explicit representation is well known and has 
the form 

X y z 

A[G]{x,y,z)^ jGi{^,yo,zo)d^ + j G2{x,Q, zo)dC + j G^ {x,y,r])dri + C. 

xo Vo zo 

Consider the equation 

(-A + I/) 5 = in G (46) 

q2 -^2 q2 

where A = + + v and g are complex valued functions, and G is a 
domain in K^. We assume that g is twice continuously differentiable. 

Theorem 40 Let f be a nonvanishing particular solution of j^6p . Then for any 
scalar twice continuously differentiable function g the following equality holds, 

{D + M^){D~M^)g^{-A + iy)g. (47) 

Proof. This is a direct calculation based on the Leibniz rule 1)44(1 . ■ 

Remark 41 The factorization was obtained in J^, ^ in a form which 
required a solution of an associated biquaternionic Riccati equation. In 
it was shown that the solution has necessarily the form Df / f with f being a 
solution of {-i&il - 
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Remark 42 Theorem \4Ul generalizes theorem^ In a two-dimensional situation 
reduces to ^10\) . 

Remark 43 As g in j^Tj ) is a scalar function, the factorization of the Schrodinger 
operator can be written in the following form 

{D + M^)fD{f-^g) = {-^ + v)g, 

from which it is obvious that if g is a solution of {4-6^ then the vector F = 
fD(f'^^g) is a solution of the equation 

{D + M^)F =^0 mG. (48) 

The inverse result we formulate as the following statement. 

Theorem 44 Let F be a solution of o, simply connected domain G. 

Then g = fA[f^^F] is a solution of i 

Proof. First, in order to apply the operator A to the vector f~^F we should 
ascertain that indeed, 

rot(/-iF) = 0. (49) 
For this, consider the vector part of H48I) . It has the form 

rot F + [F X ^y-] = 

which is equivalent to equation H49|) . 

Now, applying the Laplacian to g — fA[f^^F] and taking into account that 
/ is a solution of H46|l and F is a solution of H48() we obtain the result: 

-A.g - D^g = D{Df ■ A[f-'F] + F) 
= f-'FDf -A[f-'F]Af + DF 

= F^-ufA[f-'F]-F^ 

= -vg- 



In the same way as in section 13 we obtain the factorization of the operator 
divpgrad+g where div and grad are already operators with respect to three 
independent variables. 

Theorem 45 Let uq be a nonvanishing particular solution of the equation 

(div p gr ad+g)ii = inGcR^ (50) 

with p, q and u being complex valued functions, p € C^{G) and p in G. 
Then for any scalar function ip d G'^ (G) the following equality holds 

(divpgrad+g)^ = -~p^^^{D + M^){D ~ M^)p^/'^ip (51) 
where f = p^^'^UQ. 
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Proof. This is analogous to the proof of theorem O ■ 
Thus, if u is a solution of equation (|50() then 

F = fDif-'p'^'u) = fD{u^\) 

is a solution of equation H48|) (see remark 031 • The inverse result has the fol- 
lowing form. 

Theorem 46 Let F be a solution of equation *^ ^ simply connected domain 
G, where f = p^^'^uo and Uq be a nonvanishing particular solution of \5U\) . Then 

u = UQA[f-^Y] 

is a solution of \5U(l . 

Proof. This is a corollary of theorem|^and relation (divpgrad +q) = p^/^(A — 
u)p^/'^ where V = Af / f. m 

Notice that due to the fact that in H51|) (p is scalar, we can rewrite the 
equality in the form 

(divpgrad+<7)95 = V^^P + M^){D - ^Ch)p^''^^. 



Now, consider the equation 



Df 

{D - -j-CH)W = 0, (52) 



where W is an H(C)-valued function. Equation H52I) is a direct generalization of 
the Vekua equation (|18|l . Moreover, we show that it preserves some important 
properties of ((TH|l . 

Theorem 47 Let W = Wq + W &e a solution of f^) . Then Wq is a solution 
of the stationary Schrodinger equation 

-AWo + vWo = 0, (53) 

where v = A///. Moreover, the function u — /~^Wo fl solution of the 
equation 

div(/2gradM) = (54) 
and the vector function v = /W is a solution of the equation 

rot(/"2rotv) = 0. (55) 

Proof. Equation H52() is equivalent to the system 

divW + <^^,W^ 0, 
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Vf 

Vl^o - -jWo = 



rotW 

which can be rewritten in the form 

div(/W) = 0, (56) 

rot(/W) + / grad(/-iTyo) = 0. (57) 

From (|57|l we obtain H54|) and (|55|) . Equation (|53|) is obtained from (I54|) and 
(O. ■ 

Remark 48 Observe that the Junctions 

i j k 

Fo = f, Fi — F2 — —, F3 = J 

give us a generating quartet for the equation \5^) . They are solutions of \5^) 
and obviously any '^.{'C) -valued function W can be represented in the form 

3 

3=0 

where ipj are complex valued functions. It is easy to verify that the function W 
is a solution of ^5^1 iff 

3 

Y,{D^,)F,=0 (58) 

in a complete analogy with the two-dimensional case (see remark\^. Denote 

W ^ ipo + (fill -\- (p2j + <P3k. 
Then 1^58\) can be written as follows 

D{w + w)f + D{w~w)j =0 

which is equivalent to the equation 

Dw = YTp^"^- 

Remark 49 The results of this section remain valid in the n-dimensional sit- 
uation if instead of quaternions the Clifford algebra CZo,n (see, e.g., J^, 
is considered. The operator D is then introduced as follows D = X]j=i6iaf~ 
where Cj are the basis elements of the Clifford algebra. 
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8 Conclusions 



We showed the possibUity of a factorization of the operator (divpgrad+(7) and 
investigated only some of its appUcations. In a two-dimensional situation under 
quite general assumptions a complete system of null solutions of the operator 
can be constructed explicitly. It is quite possible that in a multidimensional case 
using the results of the preceding section the same can be done. This requires 
a multidimensional generalization of L. Bers' theory of formal powers. 

Another open question is the proof of expansion and convergence theorems 
for the bicomplex Vekua equation of the form H18|l . 
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